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ABSTRACT 
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Nuclear and starburst activity are known to often occur concomitantly. //er.st7ie/-SPIRE provides sampling of the far-infrared (FIR) spectral energy 

distributions (SEDs) of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study 

(-H large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields 

O , observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. 

I We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5cr detections at 250 pm, in agreement with previous 

O ' (sub)mm AGN studies. Their combined Spitzer-MXPS and //eric/ieZ-SPIRE colours (specifically S25o/Sio vs Sjo/Sii) quite clearly separate them 

ri ■ from the non-AGN, star forming galaxy population, as their 24 pm flux is dominated by the hot torus emission. However, their SPIRE colours 

C/3 J alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their 

Cu FIR emission. For objects at z > 2 we find a correlation between the infrared luminosity attributed to the starburst component, Lsb, and the AGN 

accretion luminosity, Lace, with Lsb '^ ^acc^- Type 2 AGN detected at 250 //m show on average higher Lsb than type 1 objects but their number is 

, ' still too low to establish whether this trend indicates stronger star formation activity. 

^ Key words, galaxies: active - galaxies: Seyfert - galaxies: star formation - infrared: general - quasars: general 

^ ' 1. Introduction and it is only now, with the advent of SPIRE (IGriffin et al. 20T0l l. 

that the full FIR SEDs of galaxies can be built, all the way to 
l/~j Active galactic nuclei (AGN) and starburst activity, both among cqq 

O , the most energetic extragalactic phenomena, have been studied j^^ ^^^j^ presented here intends to build on the previous ex- 

O . separately for decades. However, it is only in the past decade and ^^^-^^^^ ^^^^^^^^^ ^-^^^ g^-^^^^ ^^^ BLAST. We study large sam- 

^t; ■ a half, with the advent of the i nfrared (IR) o bservatories (ISO, pj^^ ^f spectroscopically confirmed type 1 and type 2 AGN in the 

> . Spitzer and now Herschel; Pilbr att et al. 201 0D that it has become j Herschel Multi-tiered Extragalactic Survey (HerME^ 

- - clear that the two phenomena are related and, more o ften than q^-^^^ ^^ ^j -^ ^^^^^ ^^j^^ observed during the science demon- 

not, happen concomitantly (e.g.g^ hweitzer et al. 20 06). ^^^^^-^^ pl^^^^ (S^P). The SEDs and IR properties of many 

^ The Spitzer IRAC and MIPS cameras sampled the spectral ^f j^e objects in the sam ples have been studied in the past 

energy distribution (SED) of extragalactic sources exacdy where (g g Richards et al .^006l Hatziminaoglou et al. 2008, 2009) and 

the peak of the AGN dust emission is expected to occur under the j^eir AGN properties are well constrained. The idea is to extend 

paradigm of an axisymmetnc dust distribution (often referred to ^^j^h study to the larger wavelengths now observed by SPIRE 

as torus) surrounding the central super-massive black hole. The ^^ ^^ ^f^^^ to ^jso constrain the starburst component of these 

IRS spectrograph provided the necessary details of the silicate objects 
feature in emission and in absorption and we now have a bet- 
ter understanding of the physics of hot dust around AGN (e.g. 

IHao et al. 2007tlLevenson et al. 20071 . The peak of the cold dust 2. AGN samples, their SPIRE detections and 
emission, however, a tracer of star formation, was beyond the mid-tO-FIR colours 

wavelength range explored by Spitzer for high redshift sources 

The AGN master sample considered here consists of a total of 

* Herschel is an ESA space observatory with science instruments 469 spectroscopically confirmed type 1 and type 2 AGN in the 

provided by European-led Principal Investigator consortia and with im- 

portant participation from NASA ' http://hermes.sussex.ac.uk 
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Lockman-SWIRE (LS) and the Spitzer First Look Survey (FLS) 
fields, with redshifts that extend to z > 4. More specifically, it 
includes SDSS quasars in LS and FLS; mid-infrared (MIR) se- 
lected AGN in FLS with spectroscopic redshifts from the MMT- 
Hectospec (Papovich et al. 2006); two MIR-selected type 2 AGN 
samples from Martinez-Sansigre et al. (2006) and Lacy et al 
(2007); and two X-ray selected type 2 samples from Mainieri et 
al. (2002) and Polletta et al. (2006). Their SPIRE fluxes are es- 
timated via linear inversion methods, using, whenever available, 
the positions of known 24 yum sources as priors (see Roseboom et 
al. in prep.). The details of the various subsamples and the num- 
bers of 250 /^m and additional 350 yum 5cr detections are given 
in Table [1] Note that none of the X-ray selected type 2 AGN 
from Mainieri et al. (2002) were detected by SPIRE. All 5cr de- 
tections at 250yL(m and 350;um have fluxes above 12.8 mJy and 
12.2 mJy, respectively. No flux cut is applied at 500 //m, these 
fluxes are used, instead, at face value. The different detection 
rates of the various samples in the two fields are the result of the 
different depths (FLS being deeper than LS, OUver et al. 2010b, 
this volume) as well as the selections of the various samples. 



Table 1. AGN sample reference, field, type, number of objects, 
objects with 5cr detections at 250 /im and at both 250 and 350 



Sample 



Type Nobj 5cr250yum 



SDSS (LS) 

SPSS (FLS) 

Papovich et al. 2006|(FLS) 

Lacy et al. 2007 (F LS) 

Martinez-Sansigre et al. 2006| (FLS) 

glletta et al. 2006](LS) 

Mainieri et al. 2002 '(LS) 
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Fig. 1. MIPS-SPIRE colours (a) and SPIRE colours (b) for type 1 
(circles) and type 2 (large squares) AGN as a function of redshift, 
over-plotted on the entire SPIRE catalogue in the FLS field, for 
objects with 5cr detections at 250 fim (a) and 350|/m (b). 



Total 



469 



156 (96) 



It is worth noting that one third of the objects have 
5cr detections at 250 jum, also reported by |Elbaz et al. (2010)| 
This detection rate is in excellent agreement with results ob- 
tained from the studies of bright (Mb < -26.1), high red- 
shift (z > 1.8) quasars at 1.2 mm, with the MAx-Planck 
Millimeter BOlometer (MAMBO) array at IR AM 30-metre 
telescope (iCarilliet al. 2001 , Omont et al. 2001J. It is, however, 
higher than the ~I5%, reported by ,Priddey et al. (2003) for 
1.5 < z < 3.0 radio-quiet, luminous (Mb < -27.5) quasars, 
observed at 850 //m with SCUBA on the James Clerk Maxwell 
Telescope (JCMT). Even though our derived detection rate holds 
for the whole absolute magnitude range of our sample (-19 > 
Mg > -28.3) as well as in broad Mg bins (here we use Mg, as the 
SDSS g-filter is the closest to the B-band used in the other stud- 
ies), if we look at the more rare, brightest objects (Mg < -27.5) 
we only find about 20% of them having 5cr detections at 250yum, 



close to the 15% reported by Priddey et al. (2003) 



Though it is difficult to say anything about the significance 
of the detections of the rest of the objects in an objective way, 
the observed properties of the AGN with and without SPIRE 
counterparts are quite similar, both in terms of optical and IR 
fluxes and in redshift. Therefore, since the present study only 
makes use of objects with 5cr detections at 250 /im, and since 
250 yum emission arises from cold dust in the star forming re- 
gions, the 250/im-detected sample may be biased towards ob- 
jects with stronger star formation activity. The issue of SPIRE 



detection rates of AGN is addressed in more detail in Stevens et 
al. in prep. 

Figure [U shows the MIPS-SPIRE colours ^aso/^vo vs 
570/^24 (a) and SPIRE colours Ss^io/Siso vs S350/S250 (b) of 
all AGN with 250 //m and 350jt/m detections above 5cr (type 
1: circles; type 2: squares) colour-coded based on their red- 
shift, and compared to all SPIRE sources in the FLS field at 
the same SPIRE detection limits (black dots). AGN have much 
bluer colours and tend to separate nicely from the bulk of star 
forming non-active galaxies, as their 24 /vm fluxes are dominated 
by the hot dust (torus) emission. This is in agreement with re- 
cent MIR spectroscopic studies with Spitzer/IRS showing that 
the MIR quasar emission mostly arises from the torus, indepen- 
dent on their FIR properties JNetzer et al. 20071) . The presence 
of strong PAH features, however, may slightly weaken this ef- 
fect in certain redshift ranges. The 7.7 /im PAH feature, in par- 
ticular, might affect the 24 //m flux of z ~ 2 AGN, when present 
(ILutz et al. 2008l l. 

In the SPIRE bands the colours of AGN are indistinguishable 
from those of the star forming non-active galaxies, suggesting 
that AGN actually appear like starburst galaxies in the FIR. This 
is consistent with longer wavelength studies that show the FIR 
emission of submm-luminous quasars and their starburst to share 
many properties of submm galaxies. Their FIR luminosities are 
similar, and the best evidence is probably that they follow prac- 
tically the same FIR/CO luminosity relation as submm galaxies 
and local ULIRGs (see e.g. Fig. 5 of lRiechers et al. 20061 ) . 
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A variety of multi-wavelength datasets covering the en- 
ergy range from the radio to the X-rays is available for 
the objects in our master sample, but in order to have a 
uniform wavelength SED coverage to perform the fitting, 
we only used the available photometry from SDSS DR7 
( |Abazajian et al. 2009|l, 2M ASS and 2MASSx6 in the Lockman 
lield (Beichman et al. 2003'), Spitzer IRAC and MIPS data from 
the SWIRE (Lonsdale et al. 2004) an d ELS Spitzer Survey s, and 
new SPIRE HerMES data. PACS ( |Poghtschet al. 2010| i data 
have not been used, as they are only available for parts of the 
HerMES SDP fields. 



3. Star formation in AGN 

The observed SEDs, described in Sect. |2] were compared, 
by means of SED fitting with a standard x^ minimisation, to 
a series of models comprising three diff'erent components: a 
stellar component composed by various simple stellar popula- 
tion (SSP) models build using the Padova evolutionary tracks 
JBertelli et al. 1994 ): a grid of AGN/torus models that include 
both a toroidal and a flared disk dust geometry presented in 
[Fritz et al. (2006)[ and six empirical starburst SEDs. For a full 
description of the SED fitting and individual model components 
see Hatziminaoglou et al. (2008, 2009). For reasons explained in 
detail in Hatziminaoglou et al. (2008) the reduced x^ can reach 
high values without undermining our confidence in the fits. We 



will restrict the present study to the objects with fits having re- 
duced ;t'^ < 10. 

This leaves a total of 68 (42) type 1 and 11 (7) type 2 AGN, 
with 5cr detections at 250 /im (and 350 /im), respectively. As al- 
ready mentioned, the 500 /im are taken at face value, even though 
the detection level of about one third of them falls below 2cr. The 
number drops to about 15% for the objects with a 5cr detection 
at 350 yum. Despite the low significance of some of the 500 /im 
data points, they follow nicely the observed SEDs as traced by 
the other FIR points and are unlikely to affect the fit, because 
of their large photometric eiTors. SPIRE 500 //m non-detections 
are not treated as upper limits and are excluded from the fits. 
Examples fits for a type 1 and a type 2 AGN are shown in Fig.|2l 

As a first remark we point out that, in order to reproduce 
the observed SPIRE data points, a starburst template is always 
needed, even if we allow for very large (kpc-scale) tori. Large 
tori are, in any case, unphysical in the sense that they extend 
well into the host galaxy where other physical phenomena such 
as star formation may occur, and the AGN is no longer the pri- 
mary source of dust heating. The SED fitting results in the set of 
estimations of physical parameters, describing the various com- 
ponents. Here we will focus on the accretion luminosity. Lace, the 
model luminosity of the accretion disk ranging from soft X-rays 
to the optical wavelengths scaled to the observed data points, and 
the IR luminosity of the starburst component, Lsb , integrated be- 
tween restframe 8 and 1000 //m, as a direct measure of star for- 
mation. 

Figure [3] shows Lsb as a function of Lace with each point 
colour-coded according to redshift for 250 /im 5cr detections 
alone (open symbols) and additional 350 //m detections (filled 
symbols). A broad correlation of Lsb can be seen with Lace, both 
quantities, however, also scale with redshift, as seen from the 
colour-coding of the points. 

If we divide the sample in bins of redshift, the picture be- 
comes less clear. Table |2] shows the Pearson correlation coeffi- 
cients for the quantities Lsb and Lace in five redshift bins, (the 
missing value in the last column reflects the very small num- 
ber of points in the relevant bin, for which a Pearson correlation 
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Fig. 2. Example fits for a type 1 (upper panel) and a type 2 
(lower panel) AGN. In red the observed data points; in blue 
the AGN/torus model; in light green the starburst component; 
in dark green the stellar component, and in black the total model 
SED. The reported ;t'^ are reduced. 
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Fig. 3. Lsb as a function of Lace and redshift (colour-coding) for 
all the objects (type 1: circles; type 2: large squares) with 5cr 
detections at 250 //m and 350 /^m (filled symbols) and with 5cr 
detections at 250//m alone (additional open symbols). 



would be meaningless). The very low values of the coefficients 
suggest that the observed global trend may reflect the fact that 
the sample being flux limited, more distant objects will also be 
intrinsically more luminous. A lack of correlation between Lsb 
and Lace would imply that star formation activity is not influ- 
enced by the presence of an active nucleus in the centre of the 
galaxy. Nevertheless, since the derived Lace only covers at most 
two orders of magnitude at any given redshift, and since we are 
dealing with a limited number of objects, it may be that our sam- 
ple is not adequate to detect any such correlations. Considering 
the z > 2.0 bin alone we do find a 95% (90%) probability of Lsb 
and Lace correlating as Lsb °^ ^ac?' for objects with 5cr detec- 
tions at 250//m (and additionally at 350/im). This value is in sur- 
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Fig. 4. Distribution of Lsb (integrated between 8 and 1000 yum), 
as an indicator of star formation, for type 1 (black) and type 2 
(grey) AGN with 5cr detections at 250/im. 



prisingly good agreement with that found by Wang et al. (2008) 
in their study of z ~6 quasars with MAMBO. 

Table 2. Correlation coefficients r and r' between Lsb and Lace 
in bins of objects considering the subsamples with 5cr detections 
at 250 jum alone (N) and those with both 250 and 350 yum 5cr 
detections (N'), respectively. 
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z>2.0 


18 


16 


0.496 


0.468 



torus (MIR) emission, respectively, and the middle wavelengths 
(optical-to-near IR) defining whether an AGN would be classi- 
fied as type 1 (unobscured) or type 2 (obscured). With Herschel- 
SPIRE we can now also accurately probe the cold dust compo- 
nent, heated by star formation, and look for differences between 
the types of active galaxies not only in the immediate environ- 
ments of their nuclei but also in their hosts. The somewhat sur- 
prising lack of a strong correlation between the AGN intrinsic 
properties (L^cc) and the star formation activity as traced by Lsb, 
at least up to z ~ 2, possibly an effect of the narrow span of our 
sample in Lace, requires further investigation with larger AGN 
samples that also cover a wider accretion luminosity range, in- 
cluding stacking sources undetected by SPIRE. The same holds 
for the seemingly higher star formation activity in type 2 AGN. 
Even though this may not necessarily conflict with geometric 
unification schemes, it certainly opens the way for investigations 
in new directions. 
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An important issue in the study of nuclear and star forma- 
tion activities is that of the occurrence of star formation in the 
different types of AGN. Figure |4] shows the distribution of Lsb, 
as an indicator of star formation activity, for type 1 (black) and 
type 2 (grey) AGN for objects with 5cr detections at 250/im. 
Because of the difference in number of AGN per type, the his- 
togram shows the fraction of objects per type. 

Type 2 AGN detected at 250 //m tend to have on average 
a larger Lsb, trend also observed in Fig. [3] but their numbers 
are too low to obtain a more robust estimate of this trend and 
to study the possible implications. The different criteria applied 
for the selection of type 1 and type 2 AGN (see Sect. |2]i, how- 
ever, are rather unlikely to affect this result, as none of the cri- 
teria particularly favours objects with stronger star formation: 
the SDSS quasars and the X-ray selected type 2 AGN are bi- 
ased, if anything, towards AGN-dominated objects; the Lacy at 
al. (2007) type 2 AGN sample consists of objects occupying the 
AGN-dominated part of the 5'8.o/5'4.5 vs 85,^/83,6 colour dia- 
gram (Hatziminaoglou et al. 2008); finally, the MIR selected 
type 1 and type 2 AGN have been chosen among 24 /vm emit- 
ters, but as was already shown in Fig. [T] the 24 //m flux of AGN 
is dominated by the torus component. 

4. Conclusions 

Until recently, the study of type 1 and type 2 objects focused 
mainly on their X-ray-to-MIR emission, with the two extreme 
A ranges sampling the nuclear activity (X-rays) and hot dusty 
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